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Signature of a magnetic activity cycle in HD49933 observed by CoRoT 
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Solar-like stars with an external convective envelope can develop magnetic activity cycles under the interaction of con- 
vection, rotation and magnetic fields. Even in the Sun, these dynamo effects are not yet well understood and it would 
be extremely important to extend this study to different stars with different characteristics. HD49933 is a F5V, 1.2 solar 
mass star that has been observed by CoRoT for 60 days during the initial Run and 137 more days about 6 months later. 
Thus, a total of 400 days have been covered with the two observations. Assuming that the activity cycle is proportional 
to the rotation of the star (which spins 8-9 times faster than the Sun, i.e., P ro t=3A days), CoRoT observations could be 
good to perform such study. The analysis techniques employed here hve already been successfully tested on sun-as-a-star 
observations done by VIRGO and GOLF on board SoHO and ground-based MARK-I instrument. 
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Astronomers of modern era have been observing the activity 
on the surface of the Sun since Galileo and they have been 
tracking the number of sunspots since then. The 1 1-year cy- 
cle of the observed active regions on the solar surface is a 
consequence of a dynamo process running in the convective 
region of the Sun (Dikpati & Gilman 2006). However, the 
delayed onset of solar cycle 24 (e.g. Salabert et al. 2009) 
gives evidence of a lack of understanding the physical pro- 
cesses governing the activity cycle. Thus, observing such 
activity cycles in other stars would, in principle, help to im- 
prove our knowledge in this field. 

Stellar activity of hundreds of stars have already been 
monitored over the last 40 years with the Mount Wilson 
survey (Wilson 1978; Baliunas & Vaughan 1985; Baliunas 
et al. 1995) or the Lowell Observatory survey (Hall, Lock- 
wood & Skiff 2007). Many stars seem to have periodic cy- 
cles between 2.5 and 25 years. Moreover, it has been sug- 
gested that there are two different branches of stellar cy- 
cles in solar-type stars, one active and one inactive (Saar 
& Brandenburg 1999), which could reflect the existence of 
two different dynamos operating in their interiors. 

To properly study the dynamo models, we need to com- 
plement the classical picture of the stars with asteroseismic 
informations (e.g. Stello et al. 2009). This will not only give 
access to complementary information on the activity cycle 
-due to the changes in the oscillation mode frequencies and 



amplitudes- but also provide direct information of the struc- 
ture and dynamics of the stellar interior. In particular, it 
would allow us to study the extension of the convective zone 
(e.g. Ballot, Turck-Chieze & Garcia 2004), the differential 
rotation, and distinguish between surface rotation rate and 
the stellar inclination axis whenever possible (see the dis- 
cussion on Gizon & Solanki (2003) and Ballot, Garcia & 
Lambert (2006)). 

In the Sun, the central frequency of the acoustic modes 
(p modes) changes between the maximum and the minimum 
of the solar cycle: the modes are shifted by «0.45 /iHz to- 
wards higher frequencies at the maximum of the magnetic 
activity (e.g. Chaplin et al. 2007). Indeed, not only the fre- 
quencies, but all the characteristics of the p modes change 
with the activity cycle (e.g. Jimenez-Reyes, et al. 2007). 
In particular, these frequency shifts -during a maximum of 
activity- are also accompanied in the Sun by a reduction of 
the observed amplitude of the modes in both: velocity and 
intensity fluctuations (e.g. Jimnez-Reyes et al. 2003). 

For a cool star, like the Sun, with an alpha-Omega dy- 
namo, a longer period rotation implies a longer cycle period 
(e.g. Thomas & Weiss 2008): 
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with R =P ro t i V c , the Rosby number and r c the convective 
turnover time and q changing from 0.25 to 1 (e.g. Ossendri- 
jver 1997; Saar 2002; Jouve, Brown & Brun 2010). Using 
t c « 12d and following the computations done by Kim 
& Demarque (1996), we find that R /R q w 0.4, hence 
P cyc » 200 - 400 days. 
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The NASA Kepler mission (Koch et al. 2010) is going 
to provide a unique opportunity to do asteroseismic mea- 
surements of the same field during more than 3.5 years of 
solar-like stars (Chaplin et al 2010), red giants (Bedding 
et al. 2010), stars belonging to open clusters (Stello et al. 
2010), and even discover new hybrid pulsators (Grigahcene 
et al. 2010). In all cases we will be able to study long pe- 
riod activity cycles (e.g. Karoff et al. 2009) as well as their 
seismic properties. This was impossible from ground due 
to the shortness of the observational campaigns (e.g. Ar- 
entoft et al. 2008). We already have more than 3 years of 
data coming from the CoRoT mission (e.g. Michel et al. 
2008) in which several solar-like stars have been observed 
(Barban et al. 2009; Garcia et al. 2009; Mosser et al. 2009; 
Deheuvels et al. 2010; Mathur et al. 2010a). Unfortunately, 
the longest observations available are limited to around 150 
days. There is one exception, the solar-like star HD49933. 
Indeed, assuming that magnetic activity cycles in stars are 
proportional to its rotation, HD49933 could be a good target 
to seismically study the existence of a stellar activity cycle. 
Its surface spins around 8-9 times faster than the Sun (Ap- 
pourchaux et al. 2008; Benomar et al. 2009) which means 
that the time span of 400 days covered by the two observa- 
tion runs could cover a significant part of such a magnetic 
cycle (Garcia et al. 2010). This star was indeed known to be 
active (Mosser et al. 2004). 

2 Methodology and data analysis 

To check whether or not there is a magnetic activity cycle in 
HD49933 we have performed three different, but comple- 
mentary studies: a) study the evolution of the surface rota- 
tion peak in the PSD, b) look for a frequency shift in the p 
modes, and c) study the rms maximum amplitude per radial 
mode, A max . 

We have analyzed the so-called '"Helreg" (Auvergne et 
al. 2009) light curves of HD49933 provided by the CoRoT 
satellite during ^60 days of the initial run in 2007 (Appour- 
chaux et al. 2008) and during ~137 more days in 2008 sepa- 
rated about 6 months (Benomar et al. 2009), covering a total 
of 400 days. 

3 Time evolution of the active regions 

The passage of the active regions crossing the visible stellar 
disk is a very good indicator of the activity. Mathur et al. 
(2008) used this methodology to study the activity cycle of 
the Sun using helioseismic radial velocity observations of 
the GOLF/SoHO instrument (Gabriel et al. 1995; Garcia et 
al. 2005) compared to 5 more classical activity indexes as 
well as to see the seismic signature after solar flares (Ku- 
mar et al. 2010). Apart from a small time lag between 3 to 
5 days, the correlation with the classical activity index was 
above 50% reaching in one of the cases ^60%. Thus, we 
start by inspecting the light curve of the combined obser- 
vations in which the presence of starspots is clearly visible 



(see Top panel of Fig. 1). We use a Global Wavelet Power 
Spectrum (GWPS, Torrence & Compo 1998) to study the 
energy distribution with time of the surface rotation peak 
at low frequency corresponding to a period of 3.4 days as 
measured by Appourchaux et al. (2008). 

We first confirm the surface rotation rate of ~3.5d with 
some differential rotation between 2 and 5 days (see middle 
panels in Fig. 1). If we divide the observations in three re- 
gions of 60 days -denoted by A, B, and C- the GWPS shows 
a quiet period -less perturbations due to crossing starspots- 
in B while there are much more GWPS power at the be- 
ginning of A. This analysis is confirmed by collapsing the 
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Fig. 1 Top: Light curve of HD49933. Middle: GWPS and 
collapsed spectrum in the time direction showing a main 
surface rotation rate of ^3.5 days. Bottom: Collapsed spec- 
trum in the frequency direction. 



4 Study of the frequency shift 

When the precision in the determination of the frequencies 
is not high enough, the frequency shift of the p modes can 
be obtained by analyzing the cross-correlation function be- 
tween the power spectra of different periods of observa- 
tion centered on the p-mode hump (Palle, Regulo & Roca 
Cortes, 1989). To do so, three power spectra of 60 days (res- 
olution of 0. 19 //Hz) have been calculated using the two ob- 
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servational runs of HD49933. To reduce the noise, the con- 
vective background modeled with a standard Harvey model 
has been subtracted (e.g. Lefebvre et al. 2008). This back- 
ground has three components and six free parameters (ig- 
noring the p-modes hump). The cross-correlation function 
is thus computed in the region of the p-mode hump between 
1000 and 2500 //Hz. 

Cross-correlation functions (see Fig.|2]) look like Gaus- 
sians. To determine a possible shift in this cross-correlation 
function, we use a frequency range of 7/iHz around the 
maximum of the cross correlation to estimate the third order 
moment of this function, which measures the asymmetry. 
Then, starting with the lag given by the asymmetry, we fit a 
Gaussian function on a window of ±5/xHz, that is approx- 
imately the second order moment of the cross-correlation 
function. The maximum of this fitted Gaussian is used as 
the position of the cross-correlation peak. In the case of so- 
lar data, if we use the power spectrum of the minimum as a 
reference, a positive shift is found. 

The result of applying this technique to the three inde- 
pendent segments of 60 days is shown in Fig. [2] The cross- 
correlating between B and C (the two segments of the sec- 
ond run), gives a maximum of the fitted Gaussian which is 
shifted 0.17 ± 0.18 //Hz. A positive shift, although it is in- 
side the error bar. The cross-correlation of A with the first 
60 days of the long series, B, produces a shift of 0.59 ± 0.23 
fiHz. Finally, the cross-correlation of A with the second part 
of the long series, C, shows a shift of 0.23 ± 0.16 //Hz. 

Although the computed correlations are quite noisy, and 
following the lessons learnt with the Sun, they could be 
telling us that a maximum activity occurs during the first 
short run A, then a minimum arrived during B, and there is 
an intermediate activity period during C. The same analysis 
using independent subseries of 30 days gave the same qual- 
itative results but with higher error bars due to the reduction 
of the time series. 



5 Time evolution of A max 

Finally, we use the A2Z pipeline (Mathur et al. 2010b) to 
obtain the rms maximum amplitude per radial mode, A max , 
of the p-mode hump and its variation with time. To do so, we 
cut the time series in subseries of 20-day long shifted by 5 
days. The convective background is fitted in each subseries 
with a Harvey-law function as it has been already described 
to compute the frequency shifts. A max is computed by fit- 
ting a Gaussian function in the frequency region between 
1000 and 2300 //Hz. 

We found a maximum in A max (see Fig. [3j around the 
middle of the long run (between B and C). This would in- 
dicate a minimum activity, which is in agreement with the 
frequency shift found when using cross-correlations. How- 
ever, it is a little shifted in time compared to the quiet period 
observed in the light curve. 



5.0 r 

4.5 - 




nn f i i i 1 

100 200 300 400 

Time (days) 

Fig. 3 Time evolution of the rms maximum amplitude per 
radial mode. The horizontal error bars represent the range 
of each subseries of 20 days. 

The amplitude during the first period, A, has an inter- 
mediate amplitude while, according to the frequency shifts, 
it should closer to a minimum amplitude. 



6 Conclusions 

We found a significant shift in the average p-mode frequen- 
cies of HD49933 as well as a variation in the maximum am- 
plitude per radial mode with time. Using the Sun as a refer- 
ence, the frequency-shift found in the analysis is telling us 
that a period of minimum activity arrives during the first 
part of the second run. A visual inspection of the light- 
curves also reveals the existence of a magnetic activity cy- 
cle because the signature of the perturbations induced by 
the active regions on the stellar surface changes with time. 
However, it is not possible to obtain direct inferences on the 
activity cycle from this observable because the inclination 
angle of the star should be taken into account as well as 
some hypothesis on the active longitudes of the star. From 
the results obtained analyzing the collapsed spectrum in the 
frequency direction, and the results from the variation of the 
maximum amplitude per radial mode, this minimum in the 
activity seems to occur at the end of the first part of the long 
run. The period of high activity in the short run unveiled by 
the analysis of the frequency-shift is confirmed by the col- 
lapsed global wavelet spectrum, and from the result of the 
variation of A max . Moreover, the frequency shifts and the 
variation of ^4 max are antic orrelated as it is the case in the 
Sun. All the analyzed information indicates the existence of 
a magnetic activity cycle in HD49933. 
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Fig. 2 Cross-Correlations functions computed between the three 60-days segments A, B and C as described in the text. 
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